Overdestabilization vs Overstabilization in the Theoretical Analysis of f-Orbital Covalency.

The complex nature of the f-orbital electronic structures and their interaction with the chemical environment pose significant computational challenges. Advanced computational techniques that variationally include scalar relativities and spin-orbit coupling directly at the molecular orbital level have been developed to address this complexity. Among these, variational relativistic multiconfigurational multireference methods stand out for their high accuracy and systematic improvement in studies of f-block complexes. Additionally, these advanced methods offer the potential for calibrating low-scaling electronic structure methods such as density functional theory. However, studies on the Cl K-edge X-ray absorption spectra of the [Ce(III)Cl6]3- and [Ce(IV)Cl6]2- complexes show that time-dependent density functional theory with approximate exchange-correlation kernels can lead to inaccuracies, resulting in an overstabilization of 4f orbitals and incorrect assessments of covalency. In contrast, approaches utilizing small active space wave function methods may understate the stability of these orbitals. The results herein demonstrate the need for large active space, multireference, and variational relativistic methods in studying f-block complexes.

What is the nature of the uranium(iii)-arene bond?

Complexes of the form [U(η6-arene)(BH4)3] where arene = C6H6; C6H5Me; C6H3-1,3,5-R3 (R = Et, iPr, tBu, Ph); C6Me6; and triphenylene (C6H4)3 were investigated towards an understanding of the nature of the uranium-arene interaction. Density functional theory (DFT) shows the interaction energy reflects the interplay between higher energy electron rich π-systems which drive electrostatic contributions, and lower energy electron poor π-systems which give rise to larger orbital contributions. The interaction is weak in all cases, which is consistent with the picture that emerges from a topological analysis of the electron density where metrics indicative of covalency show limited dependence on the nature of the ligand - the interaction is predominantly electrostatic in nature. Complete active space natural orbital analyses reveal low occupancy U-arene π-bonding interactions dominate in all cases, while δ-bonding interactions are only found with high-symmetry and electron-rich C6Me6. Finally, both DFT and multireference calculations on a reduced, formally U(ii), congener, [U(C6Me6)(BH4)3]-, suggests the electronic structure (S = 1 or 2), and hence metal oxidation state, of such a species cannot be deduced from structural features such as arene distortion alone. We show that arene geometry strongly depends on the spin-state of the complex, but that in both spin-states the complex is best described as U(iii) with an arene-centred radical.

Volatile lanthanide complexes with fluorinated heptadentate ligands.

Understanding factors that influence the volatility of lanthanide complexes remains an important goal for applications such as gas-phase f-metal separations and the synthesis of lanthanide-containing thin films. Lanthanide complexes often exhibit volatility differences that depend on the ability of ligands to saturate the lanthanide coordination sphere and attenuate intermolecular bonding in the solid state. This can make it difficult to assess how electronic factors associated with differing ligand substituents influence volatility. Here we describe the synthesis, structures, and thermal properties of a series of volatile lanthanide complexes (Ln = Nd, Er, and Yb) containing N4O3 ligands decorated with different alkyl and fluoroalkyl substituents (CF3, CF2CF2CF3, Me, and tBu). These ligands completely enveloped the tested lanthanides to form monomeric complexes with 7-coordinate distorted capped octahedral coordination geometries, as determined using single-crystal X-ray diffraction. Thermogravimetric analysis and bulk sublimation studies show how metal encapsulation affords complexes with the same volatility regardless of metal size, even with lanthanide ions with significantly different radii such as Nd3+ and Yb3+. Most notably, the results show that increasing ligand fluorination, a strategy often used to increase the volatility of metal complexes, is not always beneficial and can significantly attenuate the volatility of lanthanide complexes depending on location with respect to other substituents in the ligand framework. A pair-wise model based on density functional theory shows that the net intermolecular interactions in the unit cell can still be stronger when fluorination is present. In other words, even if individual interactions between neighboring molecules are weaker, the total number of interactions in the solid arising from the nature of crystal packing is equally important to consider.

The Influence of Phosphorus Substituents on the Structures and Solution Speciation of Trivalent Uranium and Lanthanide Phosphinodiboranates.

Here, we report the mechanochemical synthesis and characterization of homoleptic uranium and lanthanide phosphinodiboranates with isopropyl and ethyl substituents attached to phosphorus. M(H3BPiPr2BH3)3 complexes with M = U, Nd, Sm, Tb, and Er were prepared by ball milling UI3(THF)4, SmBr3, or MI3 with three equivalents of K(H3BPiPr2BH3). M(H3BPEt2BH3)3 with M = U and Nd were prepared similarly using K(H3BPEt2BH3), and the complexes were purified by extraction and crystallization from Et2O or CH2Cl2. Single-crystal XRD studies revealed that all five M(H3BPiPr2BH3)3 crystallize as dimers, despite the significant differences in metal radii across the series. In contrast, Nd(H3BPEt2BH3)3 with smaller ethyl substituents crystallized as a coordination polymer. Crystals of U(H3BPEt2BH3)3 were not suitable for structural analysis, but crystals of U(H3BPMe2BH3)3 isolated in low yield by solution methods were isostructural with Nd(H3BPEt2BH3)3. 1H and 11B NMR studies in C6D6 revealed that all of the complexes form mixtures of monomer and oligomers when dissolved, and the extent of oligomerization was highly dependent on metal radius and phosphorus substituent size. A comprehensive analysis of all structurally characterized uranium and lanthanide phosphinodiboranate complexes reported to date, including those with larger Ph and tBu substituents, revealed that the degree of oligomerization in solution can be correlated to differences in B-P-B angles obtained from single-crystal XRD studies. Density functional theory calculations, which included structural optimizations in combination with conformational searches using tight binding methods, replicated the general experimental trends and revealed free energy differences that account for the different solution and solid-state structures. Collectively, these results reveal how steric changes to phosphorus substituents significantly removed from metal coordination sites can have a significant influence on solution speciation, deoligomerization energies, and the solid-state structure of homoleptic phosphinodiboranate complexes containing trivalent f-metals.

A Perspective on Sustainable Computational Chemistry Software Development and Integration.

The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.

Molecular Geometry and Electronic Structure of Copper Corroles.

Copper corroles are known for their unique multiconfigurational electronic structures in the ground state, which arise from the transfer of electrons from the π orbitals of the corrole to the d-orbital of copper. While density functional theory (DFT) provides reasonably good molecular geometries, the determination of the ground spin state and the associated energetics is heavily influenced by functional choice, particularly the percentage of the Hartree-Fock exchange. Using extended multireference perturbation theory methods (XMS-CASPT2), the functional choice can be assessed. The molecular geometries and electronic structures of both the unsubstituted and the meso-triphenyl copper corroles were investigated. A minimal active space was employed for structural characterization, while larger active spaces are required to examine the electronic structure. The XMS-CASPT2 investigations conclusively identify the ground electronic state as a multiconfigurational singlet (S0) with three dominant electronic configurations in its lowest energy and characteristic saddled structure. In contrast, the planar geometry corresponds to the triplet state (T0), which is approximately 5 kcal/mol higher in energy compared to the S0 state for both the bare and substituted copper corroles. Notably, the planarity of the T0 geometry is reduced in the substituted corrole compared with that in the unsubstituted one. By analyzing the potential energy surface (PES) between the S0 and T0 geometries using XMS-CASPT2, the multiconfigurational electronic structure is shown to transition toward a single electron configuration as the saddling angle decreases (i.e., as one approaches the planar geometry). Despite the ability of the functionals to reproduce the minimum energy structures, only the TPSSh-D3 PES is reasonably close to the XMS-CASPT2 surface. Significant deviations along the PES are observed with other functionals.

Importance of Dispersion in the Molecular Geometries of Mn(III) Spin-Crossover Complexes.

The computational investigation of the molecular geometries of a pair of manganese(III) spin-crossover complexes is reported. For the geometry of the quintet high-spin state, density functionals significantly overestimate Mn-Namine bond distances, although the geometry for the triplet intermediate-spin state is well described. Comparisons with several wave function-based methods demonstrate that this error is due to the limited ability of commonly used density functionals to recover dispersion beyond a certain extent. Among the methods employed for geometry optimization, restricted open-shell Møller-Plesset perturbation theory (MP2) appropriately describes the high-spin geometry but results in a slightly shorter Mn-O distance in both spin states. On the other hand, extended multistate complete active space second-order perturbation theory (XMS-CASPT2) provides a good description of the geometry for the intermediate-spin state but also sufficiently recovers dispersion, performing well for the high-spin state. Despite the fact that the electronic structure of both spin states is dominated by one-electron configuration, XMS-CASPT2 offers a balanced approach, leading to molecular geometries with much better agreement with experiment than MP2 and DFT. A scan along the Mn-Namine bond demonstrates that for these complexes coupled cluster methods (i.e., DLPNO-CCSD(T)) also yield bond distances in agreement with experiment while multiconfiguration pair density functional theory (MC-PDFT) is unable to recover dispersion well enough, analogous to single-reference DFT.

Multifunctional Amine Modifiers for Selective Dehydration of Methyl Lactate to Acrylates.

Dehydration of methyl lactate to acrylic acid and methyl acrylate was experimentally evaluated over a Na-FAU zeolite catalyst impregnated with multifunctional diamines. 1,2-Bis(4-pyridyl)ethane (12BPE) and 4,4'-trimethylenedipyridine (44TMDP), at a nominal loading of 40 wt % or two molecules per Na-FAU supercage, afforded a dehydration selectivity of 96 ± 3% over 2000 min time on stream. Although 12BPE and 44TMDP have van der Waals diameters approximately 90% of the Na-FAU window opening diameter, both flexible diamines interact with internal active sites of Na-FAU as characterized by infrared spectroscopy. During continuous reaction at 300 °C, the amine loadings in Na-FAU remained constant for 12BPE but decreased as much as 83% for 44TMDP. Tuning the weighted hourly space velocity (WHSV) from 0.9 to 0.2 h-1 afforded a yield as high as 92% at a selectivity of 96% with 44TMDP impregnated Na-FAU, resulting in the highest yield reported to date.

Samarium Diiodide Acting on Acetone-Modeling Single Electron Transfer Energetics in Solution.

Samarium diiodide is a versatile single electron transfer (SET) agent with various applications in organic chemistry. Lewis structures regularly insinuate the existence of a ketyl radical when samarium diiodide binds a carbonyl group. The study presented here investigates this electron transfer by the means of computational chemistry. All electron CASPT2 calculations with the inclusion of scalar relativistic effects predict an endotherm electron transfer from samarium diiodide to acetone. Energies calculated with the PBE0-D3(BJ) functional and a small core pseudopotential are in good agreement with CASPT2. The calculations confirm the experimentally measured increase of the samarium diiodide reduction potential through the addition of hexamethylphosphoramide also known as HMPA.

Quantifying the Influence of Covalent Metal-Ligand Bonding on Differing Reactivity of Trivalent Uranium and Lanthanide Complexes.

Qualitative differences in the reactivity of trivalent lanthanide and actinide complexes have long been attributed to differences in covalent metal-ligand bonding, but there are few examples where thermodynamic aspects of this relationship have been quantified, especially with U3+ and in the absence of competing variables. Here we report a series of dimeric phosphinodiboranate complexes with trivalent f-metals that show how shorter-than-expected U-B distances indicative of increased covalency give rise to measurable differences in solution deoligomerization reactivity when compared to isostructural complexes with similarly sized lanthanides. These results, which are in excellent agreement with supporting DFT and QTAIM calculations, afford rare experimental evidence concerning the measured effect of variations in metal-ligand covalency on the reactivity of trivalent uranium and lanthanide complexes.

Isolation of uranium(III) primary phosphido complexes.

Low-valent uranium(III) primary phosphido complexes supported by hydrotris(3,5-dimethylpyrazolyl)borate (Tp*) were synthesized with phosphines of varying steric and electronic profiles. Compounds were characterized by multinuclear NMR spectroscopy (1H, 11B, 31P NMR), infrared spectroscopy, electronic absorption spectroscopy, X-ray crystallography, and quantum chemical calculations.

Actinide arene-metalates: 2. A neutral uranium bis(anthracenide) sandwich complex and elucidation of its electronic structure.

An unprecedented sandwich complex of the actinides is synthesized from the treatment of [UI2(HMPA)4]I (HMPA = OP(NMe2)3) (2) with 3 equiv. of K(C14H10) to give the neutral, bis(arenide) species U(η6-C14H10)(η4-C14H10)(HMPA)2 (1). Solid-state X-ray, SQUID magnetometry, and XANES analyses are consistent with tetravalent uranium supported by [C14H10]2- ligands. In one case, treatment of 1 with an equiv. of AgOTf led to the isolation of U(η6-C14H10)2(HMPA)(THF) (3), formed from ring migration and haptotropic rearrangement. Complete active space (CASSCF) calculations indicate the U-C bonding to solely consist of π-interactions, presenting a unique electronic structure distinct from classic actinide sandwich compounds.

Defining the Macromolecules of Tomorrow through Synergistic Sustainable Polymer Research.

Transforming how plastics are made, unmade, and remade through innovative research and diverse partnerships that together foster environmental stewardship is critically important to a sustainable future. Designing, preparing, and implementing polymers derived from renewable resources for a wide range of advanced applications that promote future economic development, energy efficiency, and environmental sustainability are all central to these efforts. In this Chemical Reviews contribution, we take a comprehensive, integrated approach to summarize important and impactful contributions to this broad research arena. The Review highlights signature accomplishments across a broad research portfolio and is organized into four wide-ranging research themes that address the topic in a comprehensive manner: Feedstocks, Polymerization Processes and Techniques, Intended Use, and End of Use. We emphasize those successes that benefitted from collaborative engagements across disciplinary lines.

CASPT2 molecular geometries of Fe(II) spin-crossover complexes.

Using fully internally contracted (FIC)-CASPT2 analytical gradients, geometry optimizations of spin-crossover complexes are reported. This approach is tested on a series of Fe(II) complexes with different sizes, ranging from 13 to 61 atoms. A combination of active space and basis set choices are employed to investigate their role in determining reliable molecular geometries. The reported strategy demonstrates that a wave function-based level of theory can be used to optimize the geometries of metal complexes in reasonable times and enables one to treat the molecular geometry and electronic structure of the complexes using the same level of theory. For a series of smaller Fe(II) SCO complexes, strong field ligands in the LS state result in geometries with the largest differences between DFT and CASPT2; however, good agreement overall is observed between DFT and CASPT2. For the larger complexes, moderate sized basis sets yield geometries that compare well with DFT and available experimental data. We recommend using the (10e,12o) active space since convergence to a minimum structure was more efficient than with truncated active spaces despite having similar Fe-ligand bond distances.

Computational Spectroscopy of the Cr-Cr Bond in Coordination Complexes.

We report the accurate computational vibrational analysis of the Cr-Cr bond in dichromium complexes using second-order multireference complete active space methods (CASPT2), allowing direct comparison with experimental spectroscopic data both to facilitate interpreting the low-energy region of the spectra and to provide insights into the nature of the bonds themselves. Recent technological development by the authors has realized such computation for the first time. Accurate simulation of the vibrational structure of these compounds has been hampered by their notorious multiconfigurational electronic structure that yields bond distances that do not correlate with bond order. Some measured Cr-Cr vibrational stretching modes, ν(Cr2), have suggested weaker bonding, even for so-called ultrashort Cr-Cr bonds, while others are in line with the bond distance. Here, we optimize geometries and compute ν(Cr2) with CASPT2 for three well-characterized complexes, Cr2(O2CCH3)4(H2O)2, Cr2(mhp)4, and Cr2(dmp)4. We obtain CASPT2 harmonic ν(Cr2) modes in good agreement with experiment at 282 cm-1 for Cr2(mhp)4 and 353 cm-1 for Cr2(dmp)4, compute 50Cr and 54Cr isotope shifts, and demonstrate that the use of the so-called IPEA shift leads to improved Cr-Cr distances. Additionally, normal mode sampling was used to estimate anharmonicity along ν(Cr2), leading to an anharmonic mode of 272 cm-1 for Cr2(mhp)4 and 333 cm-1 for Cr2(dmp)4.

Actinide arene-metalates: ion pairing effects on the electronic structure of unsupported uranium-arenide sandwich complexes.

Addition of [UI2(THF)3(µ-OMe)]2·THF (2·THF) to THF solutions containing 6 equiv. of K[C14H10] generates the heteroleptic dimeric complexes [K(18-crown-6)(THF)2]2[U(η6-C14H10)(η4-C14H10)(µ-OMe)]2·4THF (118C6·4THF) and {[K(THF)3][U(η6-C14H10)(η4-C14H10)(µ-OMe)]}2 (1THF) upon crystallization of the products in THF in the presence or absence of 18-crown-6, respectively. Both 118C6·4THF and 1THF are thermally stable in the solid-state at room temperature; however, after crystallization, they become insoluble in THF or DME solutions and instead gradually decompose upon standing. X-ray diffraction analysis reveals 118C6·4THF and 1THF to be structurally similar, possessing uranium centres sandwiched between bent anthracenide ligands of mixed tetrahapto and hexahapto ligation modes. Yet, the two complexes are distinguished by the close contact potassium-arenide ion pairing that is seen in 1THF but absent in 118C6·4THF, which is observed to have a significant effect on the electronic characteristics of the two complexes. Structural analysis, SQUID magnetometry data, XANES spectral characterization, and computational analyses are generally consistent with U(iv) formal assignments for the metal centres in both 118C6·4THF and 1THF, though noticeable differences are detected between the two species. For instance, the effective magnetic moment of 1THF (3.74 µ B) is significantly lower than that of 118C6·4THF (4.40 µ B) at 300 K. Furthermore, the XANES data shows the U LIII-edge absorption energy for 1THF to be 0.9 eV higher than that of 118C6·4THF, suggestive of more oxidized metal centres in the former. Of note, CASSCF calculations on the model complex {[U(η6-C14H10)(η4-C14H10)(µ-OMe)]2}2- (1*) shows highly polarized uranium-arenide interactions defined by π-type bonds where the metal contributions are primarily comprised by the 6d-orbitals (7.3 ± 0.6%) with minor participation from the 5f-orbitals (1.5 ± 0.5%). These unique complexes provide new insights into actinide-arenide bonding interactions and show the sensitivity of the electronic structures of the uranium atoms to coordination sphere effects.

Plausible Emergence and Self Assembly of a Primitive Phospholipid from Reduced Phosphorus on the Primordial Earth.

How life arose on the primitive Earth is one of the biggest questions in science. Biomolecular emergence scenarios have proliferated in the literature but accounting for the ubiquity of oxidized (+ 5) phosphate (PO43-) in extant biochemistries has been challenging due to the dearth of phosphate and molecular oxygen on the primordial Earth. A compelling body of work suggests that exogenous schreibersite ((Fe,Ni)3P) was delivered to Earth via meteorite impacts during the Heavy Bombardment (ca. 4.1-3.8 Gya) and there converted to reduced P oxyanions (e.g., phosphite (HPO32-) and hypophosphite (H2PO2-)) and phosphonates. Inspired by this idea, we review the relevant literature to deduce a plausible reduced phospholipid analog of modern phosphatidylcholines that could have emerged in a primordial hydrothermal setting. A shallow alkaline lacustrine basin underlain by active hydrothermal fissures and meteoritic schreibersite-, clay-, and metal-enriched sediments is envisioned. The water column is laden with known and putative primordial hydrothermal reagents. Small system dimensions and thermal- and UV-driven evaporation further concentrate chemical precursors. We hypothesize that a reduced phospholipid arises from Fischer-Tropsch-type (FTT) production of a C8 alkanoic acid, which condenses with an organophosphinate (derived from schreibersite corrosion to hypophosphite with subsequent methylation/oxidation), to yield a reduced protophospholipid. This then condenses with an α-amino nitrile (derived from Strecker-type reactions) to form the polar head. Preliminary modeling results indicate that reduced phospholipids do not aggregate rapidly; however, single layer micelles are stable up to aggregates with approximately 100 molecules.

Improved Spin-State Energy Differences of Fe(II) Molecular and Crystalline Complexes via the Hubbard U-Corrected Density.

We recently showed that the DFT+U approach with a linear-response U yields adiabatic energy differences biased toward high spin [Mariano et al. J. Chem. Theory Comput. 2020, 16, 6755-6762]. Such bias is removed here by employing a density-corrected DFT approach where the PBE functional is evaluated on the Hubbard U-corrected density. The adiabatic energy differences of six Fe(II) molecular complexes computed using this approach, named PBE[U] here, are in excellent agreement with coupled cluster-corrected CASPT2 values for both weak- and strong-field ligands resulting in a mean absolute error (MAE) of 0.44 eV, smaller than that of the recently proposed Hartree-Fock density-corrected DFT (1.22 eV) and any other tested functional, including the best performer TPSSh (0.49 eV). We take advantage of the computational efficiency of this approach and compute the adiabatic energy differences of five molecular crystals using PBE[U] with periodic boundary conditions. The results show, again, an excellent agreement (MAE = 0.07 eV) with experimentally extracted values and a superior performance compared with the best performers M06-L (MAE = 0.08 eV) and TPSSh (MAE = 0.31 eV) computed on molecular fragments.

Buffered Coordination Modulation as a Means of Controlling Crystal Morphology and Molecular Diffusion in an Anisotropic Metal-Organic Framework.

Significant advances have been made in the synthesis of chemically selective environments within metal-organic frameworks, yet materials development and industrial implementation have been hindered by the inability to predictively control crystallite size and shape. One common strategy to control crystal growth is the inclusion of coordination modulators, which are molecular species designed to compete with the linker for metal coordination during synthesis. However, these modulators can simultaneously alter the pH of the reaction solution, an effect that can also significantly influence crystal morphology. Herein, noncoordinating buffers are used to independently control reaction pH during metal-organic framework synthesis, enabling direct interrogation of the role of the coordinating species on crystal growth. We demonstrate the efficacy of this strategy in the synthesis of low-dispersity single-crystals of the framework Co2(dobdc) (dobdc4-= 2,5-dioxido-1,4-benzenedicarboxylate) in a pH 7-buffered solution using cobalt(II) acetate as the metal source. Density functional theory calculations reveal that acetate competitively binds to Co during crystallization, and by using a series of cobalt(II) salts with carboxylate anions of varying coordination strength, it is possible to control crystal growth along the c-direction. Finally, we use zero length column chromatography to show that crystal morphology has a direct impact on guest diffusional path length for the industrially important hydrocarbon m-xylene. Together, these results provide molecular-level insight into the use of modulators in governing crystallite morphology and a powerful strategy for the control of molecular diffusion rates within metal-organic frameworks.